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Abstract 



The dynamics of a driven interface in a medium with random pinning forces 
is analyzed. The interface undergoes a depinning transition where the order 
parameter is the interface velocity v, which increases as f ~ (F — FcY for 
driving forces F close to its threshold value Fc- We consider a Langevin-type 
equation which is expected to be valid close to the depinning transition of an 
interface in a statistically isotropic medium. By a functional renormalization 
group scheme the critical exponents characterizing the depinning transition 
are obtained to the first order ine = 4 — D>0, where D is the interface 
dimension. The main results were published earlier [T. Nattermann et al., 
J. Phys. II France 2 (1992) 1483]. Here, we present details of the perturba- 
tive calculation and of the derivation of the functional flow equation for the 
random- force correlator. The fixed point function of the correlator has a cusp 
singularity which is related to a finite value of the threshold F^ similar to 
the mean field theory. We also present extensive numerical simulations and 
compare them with our analytical results for the critical exponents. For e = 1 
the numerical and analytical results deviate from each other by only a few 
percent. The deviations in lower dimensions e = 2,3 are larger and suggest 
that the roughness exponent is somewhat larger than the value = e/3 of an 
interface in thermal equilibrium. 
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I. INTRODUCTION 



The ordering dynamics of an Ising magnet following a quench below the Curie (or Neel) 
temperature is controlled, to a great extent, by the motion of domain walls This is 
believed to be true also for impure ferromagnets, where the impurity may either be nonmag- 
netic (thus we are speaking of dilution or random-bond disorder) or act as to differentiate 
energetically the coexisting phases (thus we are speaking of random-field disorder). When 
the disorder is weak, either in terms of impurity strength or concentration, its effect on the 
bulk phases is minor (at least well above the lower critical dimension), but the behavior 
of the interface may be dramatically altered: the interface becomes rough in order to take 
advantage of the inhomogeneous medium, and consequently (free) energy barriers appear 
between favorable interface configurations, leading to pinning phenomena and extremely 
slow activated dynamics. 

What happens if we apply an external field F in favor of one of the coexisting phases? 
Thermo dynamically one phase prevails, but this is achievable only if the interface can get out 
of its locally favorable configurations. At zero temperature, the expected behavior for the 
interface velocity f as a function of the driving force is shown in Fig. 1. When the driving 
force is weak, the interface may adjust itself locally and settles in a nearby metastable 
position. Steady interface motion sets in only above a threshold value Fc for the driving 
force. At finite temperatures, the transition from a pinned to a moving interface is smeared 
out by thermally activated motion. At low temperatures the rounding effect is however weak 
apart from a very small region around the transition. 

The v-F curve shown in Fig. 1 is actually representative of a number of other threshold 
phenomena in condensed systems. Examples include the electrical transport of charge- 
density-waves (CDWs) and driven flux motion in type-II superconductors 0. A basic 
theoretical question in all these cases is whether one can calculate the v-F curve using a 
suitable model, which to a large extent involves an estimate for Fc and a theory for the 
singular behavior around the transition. The first part has been essentially answered twenty 
years ago in the pioneering work of Larkin |^. Progress on the second part has been made 
only recently [§,0. The difficulty has to do with the fact that the depinning transition is 
a nonequilibrium critical phenomenon which involves correlated motion of all parts of the 
system. As it turns out, the key to the solution lies in properly renormalizing the random- 
force correlator that characterizes the medium. 

The present paper is an extension of an earlier one where the basic procedure and re- 
sults of a functional renormalization group (RG) treatment of the interface problem were 
reported ^j. The strategy we followed is very similar to the ones for treating equilibrium 
critical phenomena: the interface model has an upper critical dimension = 4 above which 
perturbation theory (expanding around a flat, moving interface) converges at weak disorder. 
For D = A — e < D^, it is possible to transform the perturbation theory into a RG proce- 
dure using momentum shell integration, which yields flow equations for the parameters of 
the model. From the fixed point solution of the flow equations one obtains exponents that 
characterize the depinning transition. Our main aim here is to present details of the pertur- 
bative calculation. We shall also try to clarify the relation between the interface model used 
in our analysis and those considered by others. Finally, we compare the analytical results 
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with those from our numerical simulations [llO| , pl]| which indicate some discrepancy in lower 
dimensions. 

The paper is organized as follows. In Sec. II we introduce the model and review some of 
its basic properties, in particular the considerations that lead to the identification of an upper 
critical dimension and characteristic lengths. Section III contains formulation and results 
of a perturbation theory. The RG flow equations based on the perturbative calculation are 
presented in Sec. IV, along with a discussion of fixed point solutions and exponents. In 
Sec. V we present the results of our numerical study. Conclusions and some of the open 
problems are discussed in Sec. VI. 



II. THE MODEL AND SOME BASIC PROPERTIES 

A. The model 

On a certain coarse-grained level, the motion of an interface in an isotropic, weakly 
disordered medium in the viscous limit can be described by the following equation for the 
normal velocity f„(R) at a position R, 

Xvn = 2-fH + F + ri{R) + ifi{R,t). (1) 

Here H is the mean curvature of the interface at R, F is the driving force (which is actually a 
thermodynamic potential that couples to the order parameter), ?7(R) is a quenched random 
force, 9?(R, t) is a thermal noise, and A and 7 are friction (or viscosity) and surface tension 
coefficients, respectively. Such an equation can be motivated from a Ginzburg-Landau theory 
for the dynamics of the bulk order parameter [IS,!^,!!]- Here we simply accept (|l]) as a 



plausible dynamical model for an interface in a disordered medium that is isotropic in a 
statistical sense. 

We now choose a coordinate system with the ^-axis perpendicular to the average interface 
orientation and the remaining coordinates forming the D-dimensional x-space parallel to 
the interface. Ignoring overhangs, the interface position is uniquely specified by a function 
z{yi,t). Then Eq. (|Tp becomes 

4^ = 7V ■ (g-'/'Vz) + F + r/(x, z) + ^?-^/V(x, t), (2) 

where g = 1 + (Vz)^. On sufficiently large scales parallel to the interface, the factors g on 
the right-hand-side of (0) can be set to one. The effect of the g factor on the left-hand-side 
of (^ is more subtle. When the velocity v of the interface is finite, an expansion of g in 
the height gradient gives rise to a term ^Af (V^;)^ which changes the asymptotic scaling 
behavior to that of the Kardar-Parisi-Zhang universality class \T^\. This term vanishes as 
t> ^ 0. Thus we expect an interface model defined by |T5|-0 



dz 

A— =7V2z + F + r7(x,z), (3) 

to have the same critical behavior at the zero temperature dep inning transition as (^, 
provided the gradient of the interface height tends to zero on large scales. 
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Equation (^) is the final model to be analyzed in this paper. At this point a specification 
of the random force r] is called for. For random- field disorder, one may suppose that 77 (x, z) is 
a gaussian random variable that varies over a distance of order a (e.g., the typical extension of 
an impurity). A different situation arises for random-bond disorder, where the random force 
takes a differential form 77(x, z) = — n-VR\^(R) ^ —dV{'x, z)/dz, with l^(R) being a random 
potential correlated over a distance a and n the interface norm. A model for driven CDWs 
can also be casted in the form (^, with z being a phase variable and rj{x, z) = rjo sin[2;— /9(x)], 
where the independent random variable is the preferred phase /?(x) As far as our 

analysis goes, the three different situations allow a unified description in terms of the pair 
correlator of the gaussian random forces, 

(r/(x',z'Mx' + x,^' + z)) = A||(x) A{z), (4) 

where, without loss of generality, the mean of rj is chosen to be zero. Here and elsewhere (•) 
denotes average over the disorder distribution. The correlation of r] in x-space is understood 
to extend over a distance of the coarse-graining length ay and its precise form is not important 
in our analysis. In contrast, the functional form of the random-force correlator in 2;-space, 
A{z) = A{—z), turns out to be of crucial importance. 

Figure 2 illustrates the appropriate functional forms for A{z) for the three different cases 
discussed above. The correlator for the random- field disorder is depicted in Fig. 2a), where 
A(z) is a monotonically decreasing function of z for z > and decays exponentially to 
zero beyond a distance a±. Figure 2b) shows the correlator for the random-bond disorder, 
where A(z) = ~R"{z). Here R{z) denotes the correlator of the random potential l^(x, z) in 
2;-direction, which has the same characteristics as the random force correlator shown in Fig. 
2a). For the CDW problem [Fig. 2c)], A{z) is periodic in z and in addition the integral of 
A{z) over one period is zero. 

B. The upper critical dimension 

Suppose we accept as a reasonable model for the driven interface, what can one say 
about its properties before trying to solve it? Surely enough, one can arrive at an important 
conclusion on the existence of the depinning transition below a certain critical dimension, 
and an estimate for the threshold force Fc in the case of weak disorder, using a simple 
argument due originally to Larkin and developed in the interface context by Bruinsma 
and Aepph |T^. Let us recall the argument below. 

We begin with the remark, that a completely flat interface is never pinned. Indeed, the 
total random force on a piece of a fiat interface of extension L is roughly A^/^(0)L^/^ as 
follows easily from (4). This has to be compared with the total driving force FL^, which 
will always win for L ^ 00. 

For a rough interface this estimate is only valid on scales L at which the typical interface 
distortion h is smaller than a± and hence pinning is possible for small F. However, roughness 
of the interface will be suppressed by the curvature force '~fhL~^^^, which tends to flatten 
the interface. In D > = 4: dimension the curvature force will always win on sufficiently 
large length scales. Thus, for very weak disorder, A(0) <^ 'y'^a?~'^a']_, there is no length scale 
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down to a\\ at which the interface is rough and hence there is no pinning = 0. (Pinning 
is however possible for larger disorder). 

The situation is different for D < Df.: the pinning forces cannot be balanced by the 
elastic forces for h < a± on length scales larger than the Larkin length |[7|JT6[], 

= [7^ai/A(0)]l/^ (5) 

where e = 4 — D. Beyond this length, the rough interface explores the random environment 
to acquire a finite pinning force per unit area. The force that is needed to get the interface 
out of a typical pinned configuration is thus of the order of [|l5| , [l^] 



C. The depinning transition and scaling 

The argument presented above suggests that, below the critical dimension Dc = 4, 
there is a depinning transition at a finite strength of the driving force. An observation 
due to Middleton JlSl, originally made in the CDW context, makes it plausible that f is a 



continuous function of F. Thus the transition, if it exists, should be continuous and the 
interface is expected to exhibit critical fluctuations at the transition and crossovers away 
from the transition |^. 

We discussed already, that the interface is rough. In particular, we will assume, that 
the interface is self-affine, an assumption which has to be proven afterwards. In the ther- 
modynamic limit the height correlation function is assumed to obey the following scaling 
form 

([.(x,t) -z(x',t')]^) - |x-xf^G [j^^) , (7) 

with G{y) const for ?/ — * and 

G{y) ~ y ^ oo, (8) 

where ( and z (not to be confused with the interface height) are roughness and dynamic 
exponents, respectively. 

For F > Fc there is a length (and correspondingly, time) scale, above which one can 
neglect the non-linearity of Eq. (3), which is hidden in the random force term. Going 
over to a co- moving frame z = vt + /i(x, t), where v is the mean interface velocity and 
{h{'x,t)) = 0, we consider first the case of high velocity. Then //(x, z) can be replaced by 
?7(x, ft), i.e. 7] acts as a thermal noise. In this case the exponents are z = 2 and C = fo^' 



D > 2 and ( = — for D < 2 Upon decreasing v, the non-linearity becomes more 



important. Indeed, if a small part of the interface hits an obstacle, it will first keep its 
position, i.e. /i(x, t) in the argument of ?7(x, -|- h) will decrease by —Ah from its original 
value to compensate the increase of vt. 
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The typical value of \Ah\ increases with t according to (7) and (8) as t^/^. Since C < 1 
for a well defined interface and z > 1 (as we will show below), h changes more slowly than 
vt and there is a timescale ~ t>^^/(^~^) where the vt term starts to dominate the argument 
of rj. On this time scale, the perturbation of the interface by the obstacle has spread over a 
region of linear size 

^^v~^. (9) 

For larger time or length scales the effect of the non-linearity can be neglected and the 
above mentioned exponents apply. For smaller scales however, we expect new exponents 
dominated by the non-linear terms of (3). If the velocity decreases close to the threshold as 
a power law 

(10) 

we get C^{F - Fey" with [8] 

e = u{z-C). (11) 

Another important constraint on the exponents follows from a Harris type argument 
||20| for a homogeneous transition in a disordered medium: Indeed, let us divide the volume 
which includes the interface into correlated regions, each of linear size ^ parallel and 
perpendicular to the interface, respectively. The typical fluctuation 6Fc of the threshold 
force Fc in such a region is 

6Fe ~ F./e^^+'^^/l (12) 

To have a sharp transition, the condition 6Fc <^ F — Fc for F ^ Fc has to be fulfilled. This 
gives, together with (12), 

l/u<iD + C)/2. (13) 



D. Mean field theory 

Above the upper critical dimension Dc = 4, the dimensional analysis of Sec. II. B suggests 
absence of a depinning threshold when the disorder is sufficiently weak. This conclusion can 
be verified explicitly using a perturbation theory discussed in the next section, provided 
?7(x, z) is differentiable in z. An exception is given by the case where ^^(x, z) may increase 
abruptly as z increases (i.e., the pinning force weakens in a discontinuous way), leading to 
a sudden release of an interface segment. As the phenomenon has an analogy to avalanche 
motion in lower dimensions, let us discuss it in some detail below. 

For simplicity, let us consider a mean field approximation to Eq. (^, which can be 
obtained by taking the limit of infinite-range, suitably normalized elastic interactions or 
D oo. Then, all interface elements z(x, t) are uniformly coupled to the spatially averaged 
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position z(t) instead of a coupling between nearest neighbors as in Eq. (H). Assuming that 
the mean position z{t) moves with a constant velocity v, i.e. z{t) = vt, one obtains P, p!7| , pT 



— = j[vt-z{t)]+F + r]{z), (14) 

where we have omitted the redundant x-dependence of rj and set A = 1. 

Equation (|l^) acquires a more transparent form p2| under a new variable u = z — vt — 



(F — f )/7 which satisfies (upon a translation of rj), 

it = —^u + ri{vt + u), (15) 
with the self-consistency condition 

^(^u) = -F + v. (16) 

The driving force F now disappears from the equation of motion ([15|) and can be computed 
from ([16|) as a function of v once a solution to ([15|) is found. 

The variable u(t) in Eq. (|15|) can be interpreted as the position of an overdamped particle 
in a moving potential, 

W{u,t) = ^u^ + U{vt + u). (17) 

The random part U (vt + u) = — /J^ ^" dz r](z) travels at a constant velocity v to the left or 
to the right depending on the sign of v. 

Using the overdamped particle language, one can relate the threshold force to certain 
geometrical properties of rj without an explicit solution of Eq. ([T5|) . Let us define a threshold 
force for forward motion, Fj:~^^ = limt,^o+ and a threshold force for backward motion, 
= lim„_^o- -^('^)- For & random force with a statistical inversion symmetry [i.e. r](z) — >■ 
-v{-z)], = = and Eq. (||) yields H 



lim (u) — lim (u) . (18) 



Fisher observed that, in the limit v ^ 0, the particle stays in a local minimum of the 
potential W essentially the whole time, interrupted by occasional jumps when the minimum 
it follows disappears. When the potential travels to the left (f > 0), the particle traces 
the leftmost minimum of W (apart from an initial transient). The opposite occurs when 
the potential travels to the right. Thus Fc is simply proportional to the average distance 
between the left and right local minima of the potential W. For a continuous, slow-varying 
function ri{z) with a sufficiently small amplitude, the potential W has a unique minimum all 
the time, and hence Fc = 0. However, for a discontinuous function 7]{z) (or for sufficiently 
strong disorder), the potential W has, from time to time, more than one minimum, and 
hence a nonvanishing Fc. Note that only upward jumps in r] give rise to upward cusps in 
the potential which, at times, serve to break one minimum into two. 

The different roles played by upward and downward discontinuities of rj when the velocity 
of the particle is small can be appreciated by considering a simple, "two-state" model, which 
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has been solved explicitly . In this model, the 2;-axis is divided into a uniform sequence 



of cells, each of length a. The pinning force ri{z) is constant in each cell: rj = ±770 with equal 
probability. Both upward and downward discontinuities are thus present. When the particle 
encounters a downward discontinuity, it sticks at the cell boundary until the pulling force 
[including the vt-teim in Eq. (^^] builds up to overtake the resisting force —770 in the cell 
ahead, and then enters the new cell at a velocity of the order of v. The situation is different 
at an upward discontinuity. When the particle reaches such a point, it suddenly experiences 
a pulling force by an amount 2?7o greater, so that its velocity jumps by the same amount 
[see Eq. (|T5|) ]. This relatively fast motion makes it possible for the particle to sample the 
negative parts of rj more often than the positive parts in the limit f — > 0, giving rise to 
a finite Fc. At sufficiently small tiq (i.e., no more than two minima occurs at any given 
time), W has two minima when the upward discontinuity is within a distance tiq/j from the 
origin. The distance between the two minima is given by Au = 2?7o/7- Using Eq. ([T8|) , one 
finds Fc = 770/(207). Apart from the fast motion (i.e. a sort of instability) caused by the 
sudden increase of driving at an upward discontinuity of t], the rest of the motion is regular 
and linear response theory applies. Hence for F greater than but close to Fc one expects 
V F — Fc- The mean-field velocity exponent is thus given by 6'mf = 1- The above results 
for the threshold and for 6'mf are in agreement with the explicit solution of the model [^ 



It is interesting to see that a perturbative treatment of Eq. ( p!4D ||T7|| also yields the mean- 
field exponent 6'mf = 1, although, for the "two-state" model, it yields a threshold which is 
off by a factor of 2. In the lowest order perturbation theory, the disorder is completely 
specified by its correlator A(z). The result for the threshold is given by, 

F,^-A'(0+)/7 + O(A^), (19) 

i.e., for weak disorder, a finite threshold requires a cusp singularity of the random-force 
correlator at distance 0^. For the two-state model, A(z) = ?7q(1 — |2;|/a) for \z\ < a 
and A(z) = 0. Both upward and downward discontinuities of rj contribute to a nonzero 
A'(0+) = -vi/a. 

What can we learn from the above discussion for the case D < Dc, for which the scaling 
arguments suggest a finite threshold? We have seen that for v ^ the particle z{t) in mean 
field theory follows most of the time the minimum of its local potential. The same is true 
for an interface element z{'K,t) for D < Dc- The long periods of slow motion of z{'x,t) are 
interrupted when the local minimum disappears. This happens when elastic forces pull the 
interface element over a barrier of negative r/'s. Then the interface element rapidly moves 
forward and it can itself pull neighboring interface elements over their barriers ( "avalanche" ) . 
In mean field theory the analogous scenario is when the particle z{t) is dragged by the mean 
position z{t) = vt over an upward discontinuity of the random force. The following rapid 
motion of the particle ("jump") results in a finite threshold Fc > 0. The jumps and the 
finite threshold are a consequence of the discontinuity in ri{z), which correspond to a cusp of 
the correlator (see Eq. (|T9|)). Given the observation that the jumps (and a finite threshold) 
are present for D < Dc, one can expect that this jerky motion has to be described by an 
effective singular correlator A(z). In other words, the avalanche-like behavior for D < Dc 
can be a mechanism that gives rise to an effective discontinuous random force acting on the 
interface. We will see below in the functional RG treatment, that the fixed point function 
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of A{z) has indeed a cusp singularity at the origin, and that this cusp is responsible for a 
finite threshold, as in mean field theory p3[ ! 



III. PERTURBATION THEORY 
A. Iteration scheme 

In this section we present a perturbation theory which forms the basis of our renormal- 
ization group calculation. 

In the moving phase the interface fluctuates around its average position vt, z(x, t) = 
vt + h{-x,t). Then, Eq. (0) can be written as 

dh 

X— = ^V^/i + F-Xv + r/(x, vt + /i(x, t)). (20) 

A perturbative solution to Eq. (^) can be found by expanding 77 (x, z) around the fiat co- 
moving reference interface vt. To do so we Fourier-transform ri{:K,vt + h{'x,t)) with respect 
to its second argument and expand the term e^'^^, 

ri{x, vt + /i(x, t)) = j r7(x, q) e^'?[^*+''("'*)l 

Jq 

= J^v{^,q)e''^^\l + tqh{^,t)+ l(^g)^[Mx, t)f + ...) , (21) 

where J^^dq/{2'7i) is abbreviated by J^. Since the strength of the disorder A(0) determines 
the amplitude of the height fluctuations h, one expects that the expansion in powers of h 
converges for sufficiently weak disorder. 

Using (pID, we can rewrite (pO|) in Fourier space, 

h{k, u) = Go(k, iu){5{k)5iiu) [F - Xv] + -fj fk, -) + (^(k, 00) 
^ V \ vj 

+ 11 iq fi{k-p,q) h{p,LJ -qv) 

Jq Jp 

{tqy fj{k-p-p',q)h{p,n)h{p',uj-n-qv) + ... }. (22) 
Here, Jp = J d^p/{27T)^ and the bare propagator is defined by 



^ 2 Jq Jp Jn Jp' 



We added a source term (p(k, u) on the rhs of Eq. (|22D . The Fourier-transformed quantities, 
e.g., 

h{k,u) = J d^xdt h{^,t)e-'^^-''+'^'^ , (24) 
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are marked by a tilde. Eq. (p^) can be represented diagrammatically, as shown in Fig. 3. 
The short double arrow represents a line with an arrow is the bare propagator Golk? a 
cross X stand for 5{}<i)5{uj) [F — \v] and the open circle denotes the source term (p. A filled 
circle with a dotted line indicates the quenched noise 77, where the number of outgoing arrows 
of such a vertex determines the power of the term iq and the number oi D + 1 dimensional 
integrations. 

The averaging over the quenched noise will be represented by the connection of two 
dotted lines and is performed using 

(r/(k, g)77(k', q')) = Ay (k) A, 5^(k + k') 6{q + g'), (25) 
where A||(k) ~ e is the Fourier transform of A||(x) and Ag is defined by 

A{z) = f Ag e'^'. (26) 
Assuming A(z) is analytic at 2; = 0, we may write 

AW^Et;^^-ft„. (27) 



where 



Q2n = / g'" A, (28) 



will be called moments of the noise correlator A. 

Equation can now be solved by iteration. As an example, we insert the second 
term GoV'^fj on the rhs of Eq. (^) for the h of the fourth term which gives a "one-loop" 
contribution to the solution h{h,u). Neglecting higher order terms and averaging Eq. (^2]) 
over the disorder, we obtain (see Fig. 4) 



A A JqJp 7p^ — zAvq 



For D > Dc = 4:, the integral in the correction term converges for all v (see Appendix A.l). 

When D < 4, the integral in (p9| ) diverges in the limit v —>■ due to contributions from 
small momenta, p (see Appendix A.l). This statement can be confirmed by continuing 
the iteration to higher orders (see Appendix A. 2). Quite generally, the expansion parameter 
can be written in the form (^o/-^c)^~^5 where 

eo = iiaJAvy/' (30) 

plays the role of a bare correlation length. This parameter is small at high velocities. 
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B. Linear response and renormalization of 7 and A 



In the moving phase, the interface is free on large length and time scales. The dynamics 
of the system for small k can thus be described by linear response theory with effective 
parameters. As usual, one can define an effective linear response function G(k, u) 



dh(k, uj) 



(31) 



tfi=0 



Due to the nonlinear terms in Eq. (^21) , G(k,uj) differs from Go(k, cj). The correction can 
be expressed in terms of a "self-energy" S(k, cj). 



^(k, to) = Goik, uj) + Go(k, cj)S(k, cu)G'(k, u). 



(32) 



As we shall show, S(k, u) has the same functional form as Gg ^(k, u) at small k and u. Hence 
the correction term in (|3^) can be absorbed into the free one with the bare parameters 7 
and A replaced by effective parameters Ag// = A + 5A and 7e// = 7 + S'j. 

To calculate ^7 and 6X one has to consider diagrams with one incoming and one outgoing 
line (see Eq. (^)). The lowest order contribution to the self energy E(k, cj) is given by the 
second and third diagram on the rhs of the diagrammatic perturbation expansion for G(k, uj) 
shown in Fig. 5. The corresponding analytic expression is 



S(k,cu) - 
For small k and uj Eq. ( 
E(k, UJ 



A. 



{iqf 



q Jp 



7P^ + iAgf 7(p + k)2 + iA(u; + gi;) 
3f) can be written as 



7k2 



D 



- 1 



iXuj 



p (7P^ + iXvqY' 



(33) 



(34) 



As in Eq. (pQ]), the integration over p converges at large p as long as4 — D = e>0, while 



Xvq/'j serves as a lower cutoff for p. The q integration is truncated for large q by the 
correlator Aq which is assumed to decay as e~'^ . The integrals can be evaluated to give 



S(k,cu) ~Q2 



72 g 



1 ija^y/^ 



7 — k — iXuj 
D 



(35) 



(see Appendix A.l, where similar integrals are carried out). When comparing this with Eq. 
( p3| ) and using (^), the first order corrections 57^^-' and 5A^^^ to 7 and A are 



^7 



(1) 



K 



-7- 



5A(i) = XCQ2-C0 = A 



K 



D 



D 

Id 



(36) 
(37) 



where c = Kd/'j'^ and LI = 7^/Q2- With these notations the small parameter of the 



perturbation theory can be written as c(52^o i^o/ Lc 
the role of a coupling constant. 



We see that Q2 = -A"(0) plays 
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The result for 5A*^^^ in (|37D can also be obtained by considering the lowest order correction 
Eq. (^) to the bare velocity vq = F/X, as done by Feigel'man In Appendix A. 2 the 
calculation of the velocity is extended to second order, with the result, 



A 



A 



1 + 



5A« /5A«^ 



A 



A 



+ 



The series obtained by perturbation theory are thus divergent for v 
used directly close to the depinning transition. 



(38) 



and cannot be 



C. Interface width 

The width of the interface (/i^(x, t)) can be calculated from Eq. (|2^) in a very similar 
way as the velocity. The lowest order diagram for the width is shown in Fig. 6. The 
corresponding analytical expression is given by 



which can be evaluated easily, giving 



A. 



q Jp 



72p4 + (Xvqy 



2 2Ko if^oY 



+ 0(6) 



(39) 



(40) 



D-2 e\Lr 



0(e) 



The interface width diverges with f — > 0. Since is a correlation length, one obtains from 
Eq. ( PPD the perturbative roughness exponent (p = e/2. 

It can be shown that the second order contribution to {h'^{'x,t)) is proportional to the 
square of the first order Eq. (|40|) but with one factor 1/e lacking: 



2K 



D 



D 




(41) 



Here we have implicitly assumed that A'(0) = which will be discussed below (Sec. V.). 



D. Renormalization of the moments of A(z) 



In the preceding section we have considered the renormalization of 7 and A. Eqs. ([36| - |37D 
show that the parameter controlling the renormalization of A is proportional to the second 
moment of the disorder correlator Q2 = Jg q^^q- In order to carry out the renormalization 
group analysis of the problem it is necessary to know the renormalization prescription of Q2- 
The usual way to get this prescription is to find a vertex function, which is proportional to Q2 
when calculating the lowest-order contribution in the disorder strength. Then, the higher- 
order contributions to this vertex function yield the desired renormalization prescription of 
Q2- It turns out that a convenient quantity to look at for this purpose is 
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Gn= f /,Go^(k,^)Go^(k',a7') 



6h{k, uj) 



6ho{pi,Qi 



6h{k',uj') 



ho=o,Lo=o 6ho{p[,Q[] 

k=pi=0 



ho=0,iu'=0 
k'=p[=0 ' 



(42) 



where hQ{p,Q) = Gq{p,Q)(p{p,Q). The subscript c indicates that only connected diagrams 
are to be considered. The perturbation expansion of (^) in powers of the disorder correlator 
is obtained by iterating h(k, u) in Eq. (|^) according to Eq. ( ^2]) and carrying out the 
average over the disorder. This perturbation expansion can be represented by means of 
diagrams. The lowest order contribution to (|4^ ) in powers of A (2;) is found by replacing all 
h on the rhs of Eq. (|2^ ) by ho- The first-order contribution to ( ^21) is given by diagram 
a) in Fig. 7 with n = 1. According to Eq. ( ^21) and Eq. ( ^2]) the analytical expression 
associated with the diagram 7a) is exactly the second moment Q2- The renormalization 
prescription of Q2 can be obtained from the one-loop diagrams in Fig. 7b)-d) by setting 
n = 1. In each of these diagrams, there is only one integration over the internal momentum 
p. This integration has the form J^CKp) and yields for all diagrams under consideration 
the same factor c^o/e. In addition to this integration, there is also a factor which involves an 
integration over the momenta qi and q2 conjugated to the argument of the disorder correlator 
A(z). This factor results in 



/ / i-qiql + HqI - 3gi'g2)A(gi)A(g2) = 3Q 



(43) 



where we have taken into account that only even moments of the disorder correlator differ 
from zero. The total contribution of the diagrams in Fig. 7 for n = 1 leads to the following 
renormalization prescription of Q2, 



Q2, 



1 



eff 



g2(l + 3Q2C-a + ...)• 



(44) 



A detailed analysis of the present problem shows that, besides the renormalization of the 
second moment, higher moments Q2n {n = 2, ...) also renormalize. The renormalization of 
Q2n can be derived by considering the following generalization of the expression (52), 



Gnn= I [ ,G,\k,Lu)G^\k',u') 

I (5"^(k,cj) 



5/io(pi, fil)(5/io(p2, ^^2) ■ ■ ■ 5/io(Pn, ^n) 
(5"/l(k', CJ') 



ho=0,w=n2=...=r2„=0 
k=pi=...=p„=0 



(45) 



ho=0,w'=C'2=...=c;=0 / 

k'=p;=...=p;=o 

The perturbation expansion of ( ^5]) is obtained in the same way as the iteration of Eq. (^2[). 
The lowest-order contribution to (^3|) in powers of A (2;) is found by replacing all h on the 
rhs of Eq. ( P^ by h^. Only the diagram with n outgoing lines and one dotted line in Fig. 
3 gives a contribution to the first-order correction of As it can be easily shown, the 
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latter is equal to the bare moment Q2n (see Eq. (pT])) and is represented by the diagram a) 
in Fig. 7. The second-order contributions to Gnn in powers of A are given by diagrams b)-d) 
in Fig. 7. The renormalization of the moments Q2n and consequently the renormalization 
of the disorder correlator A(z), can be obtained by computing Gnn up to the second order 
in A. The details of such one-loop computation of Gnn are presented in Appendix B. As a 
result, the renormalization prescription of the moments Q2n is obtained as 

n 

Q2n,eff - Q2n ^ SQS. ^ E Clji^' Q2kQ2n-2k+2 'Co, (46) 

k=l ^ 

where Cf = (l) are the binomial coefficients (see Appendix B). In particular, the correction 
to the second moment Q2 is given by Eq. (0). Eq. (|46|) has been derived for n = 1,2, .... 
However, it can be extended to n = 0. Taking into account that = we get from (^61) 
that the zeroth moment Qo does not renormalize. This conclusion is in agreement with the 
result on the interface width in the preceding section. Equation (|46|) will be used in Sec. IV 
to derive the flow equation for the correlator A(z). 

IV. RENORMALIZATION GROUP 

Approaching the critical point with t; — > 0, the perturbative corrections to the parameter 
of the model diverge as v~'^^'^, which is caused by the integration over small momenta. 
A possible way to sum up such singularities is to perform a renormalization group (RG) 
procedure, where one integrates out fluctuations in an infinitesimal momentum shell only, 
leading to differential flow equations. From the flow of the relevant parameters the critical 
behavior can be deduced. The calculation is simplifled close to the upper critical dimension 
Dc = 4. We will therefore determine the critical exponents of the depinning transition by a 
RG procedure in D = 4 — e dimensions. 



A. Three-parameter RG 

The smallest closed set of flow equations to flrst order in e is given by [H 

d\n-f/d\nL = 0, (47) 
d\nX/d\nL = cQ2L% (48) 
dQ2/d\nL = 3cQlL' (49) 

(c = Kd/'j'^)- To derive these flow equations we consider 7e//, Kff and Q2,eff as depending 
on the upper cut-off length, L = ^q. Eqs. (^)-(^) can be formally obtained by differenti- 
ating Eqs. (|36|), (|37|), and (|i^ and then setting •y^ff = 7, Kff = A, and Q2,eff = Q2- This 
is equivalent to the RG transformation, where the degrees of freedom in a momentum shell 
are successively integrated out. There, the corrections to the parameters are calculated in 
the same way as in perturbation theory - the only difference is that one integrates only over 
the momentum shell instead over all momenta. 



14 



Note that the perturbative correction (|36D for 7 is finite in the hmit e — and does not 
contribute to first order to a renormahzation of 7 The flow equation p9|) can also be 
obtained from the second order correction ( p8| ) to the velocity . 

Since there is no term in the second order correction (|41| ) of the interface width, 
provided A'(0) = 0, we conclude that A(0) does not renormalize to first order in e. 

When integrating Eq. we see that Q2 becomes infinite at L ^ Lc (if we assume that 
A(z) is analytic and monotonously decaying). It then follows from Eq. (^) that also A 
diverges at L ~ Lc. This means that beyond Lc, there is no dynamic response to any finite 
force, which is an unphysical result ! 



B. Functional RG 

Eq. ( ^61) shows that besides the renormahzation of Q2 there is also a renormahzation 
of all the other moments Q2n{f^ = 2,...). The renormahzation of the moments causes a 
renormahzation of the disorder correlator A (2;) |2^. Since all Q2n have the same upper 
critical dimension we can put together the prescriptions (^) into a functional flow equation 
for A(z) (see Appendix C): 



d A(z) _ _^^e^ 

d \a.L dz^ 



(50) 



To discuss the fixed point behavior it is convenient to consider the renormalized and 
rescaled correlator A (2;), which has the original short- distance cutoff of the bare correlator. 
To do so we perform a scale transformation x Lx, t — » L^t, and h — > L^h with Eq. 
(^). At the critical point we expect scale invariance. Since 7 is not renormalized, the 
renormalized and rescaled equation of motion, 

Bh 

\L^-'^ = ^V^h + L2-C(F - \v) + L^-^7]{L^, vLH + L^h), (51) 

should be independent of L. This is the case when A(L) ~ L^~^ and 

A(z) = const. X L'^^-'K{zL-^). (52) 



To investigate the fixed point solution A*(?/) = lim2,_^oo A(y) we insert (^) into (pO|), 
take the limit L 00, and determine the constant in Eq. (0) such that A*(0) = 1, 

(e - 2C)A*(y) + Cy^*'{y) - [I\*'{y)f - A*"(y)[A*(y) - 1] = 0. (53) 

The same equation with A*(?/) = R"*{y) was already derived for the potential-potential 
correlator R{y) in the problem of an equilibrium interface subject to a random potential |^ . 
From that treatment it is known that the fixed-point function A*{y) has a singularity 
at the origin. The ansatz 

A*{y) = l + ai\y\ + ^a2y^ + ... (54) 
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satisfies Eq. ( ^31) and by a comparison of tlie coefficients of and one gets a\ = e — 2C, 
and a2 = (e — C)/3. 

Tlie development of tliis cusp singularity can be understood as follows (see Fig. 8). 
Starting the renormalization procedure at L = Lq with a smooth correlator (A"(0) < 0), 
the curvature of A(z) at the origin becomes greater under the flow (^9[). Hence, the position 
zq of the point where the second derivative A"(z) changes its sign is shifted closer to the 
origin. A diverging curvature means Zq = 0, i.e. A"{z —>■ 0^) > 0. 

Repeating the perturbative calculations with a correlator obeying Eq. (Q) we find 
that Eqs. (|3^) , (|37| ) and (^7|) , (^Sj) remain valid to first order in e. With the interpretation 
Q2 = A"(0"^) > we overcome the above unphysical result that the inverse mobility diverges 
at L = Lf^. 

Since our earlier assumption A'(0) = failed, which yields (Pl), the conclusion of 
the nonrenormalization of A(0) is no longer valid. Instead we obtain with Eq. (|54D 
dA{0)/dlnL = — cL'^[A'(0)]^ which is in agreement with Eq. (|53D for y = when using 
Eq. (H). 

The singular term \z\ yields a reduction of the driving force, F F — Fc. To see this we 
split the correlator A in a regular and a singular part. In Fourier-space, the singular part 
behaves as A^™^ ~ —q^'^A'{z 0^). Using this in Eq. (29) for the velocity, we obtain 



D-2-f 



A'(Lo,;^-0+)Ao^-^ (55) 



which depends on the upper cutoff Aq — vt/Lq of the momentum space integration. Using 
A'(0+) = A{0)/a±_ (see Eq. (|5^)) and identifying Lq with Lc yields the threshold F^ of Eq. 
(0) , in agreement with the estimate of the scaling arguments [[7|,p!5|,p!6 . 



C. Critical exponents 

The roughness exponent ( can be obtained immediately from Eq. (^), provided A*{y) 
goes sufficiently fast to zero with y 0. Integration over y yields (e — 3^) dy A*{y) = 0. 
For random-field disorder dy A*{y) > 0, from which one gets 

C = e/3. (56) 
The other exponents can be read off from Eq. (pj). To calculate the dynamical exponent 



z we have to know the scale dependence of A, which can be obtained by integrating (^ 
using (ID and (0), 

A(L) = Ao(L/Lo)-(^-^)/^ (57) 

where Aq = A(Lo). Using Eqs. (^, (|^), (p^) and imposing scale invariance of Eq. ( |HT| ) at 

F F — Fc = V = 0, z is given by 

z = 2-{e-0/3. (58) 
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As discussed earlier, above the threshold the quenched character of the noise correlator 
changes to a thermal one on length scales L > ^ ~ ^-i/(^-C)_ Hence we have to stop the 
renormalization at ^. From F — = X{L = C,)v follows 

vr^iF-F,Y, with ^ = l-li^ (59) 
with z/= 1/(2-0 (60) 



where we have used Eqs. (0) and 

Inserting the result for random-field disorder ( = e/3 in Eqs. ( ^8|) and (|59|), ( |60l) we 
obtain the final results for the critical exponents to first order in e, 

1 2 

C ~ -e, 2; ~ 2 e, 

^3 9 

. . 1 2e , 1 3 11 

fc' ~ 1 ~ 1 e, ~ — \ e. (61) 

36-e 9 6-e2 12 ^ ^ 

The roughness exponent ( = e/3+0(e^) differs from the value C,p = e/2 from perturbation 
theory (see Sec. III. C) P,p8|, but coincides with that of an equilibrium interface in a medium 



with random-field disorder |l|. Recently, Narayan and Fisher showed that higher order 
terms in Eq. (|50| ) vanish when integrated over all z. Thus, the results for the exponents 
C and V could be exact, although non-perturbative corrections cannot be ruled out 0. In 
fact, as we will see in the next section, our simulations in < 3 yield values for C, and v 
which differ slightly from the results (pT|). 

The velocity exponent 9 is smaller than one for D < A and approaches the mean field 
value 6'mf = 1 for D = 4 (see Sec. II). The superdiffusive dynamical exponent z < 2 
is consistent with the expectation that the motion of the interface close to threshold is 



governed by the behavior of avalanches [^. The Harris condition (|13D is fulfilled as an 



equality to the first order in e. Narayan and Fisher showed that the scaling relation (pOD 
between and v is exact due to a symmetry of Eq. @). Thus, the Harris condition (|13|) can 
be written as C > (4 - D)/?>. 

Narayan and Fisher investigated also the RB case (see Sec. II) and argued that the 
results ( |6TD are also valid for RB disorder. The critical exponents for CDW's can be obtained 
by inserting the corresponding value Ccdw = into the scaling relations (|5^) , (|59|) and (|60[) 
which are independent of the form of the correlator A(z) (see Sec II.) 0. 



D. Phase diagram 

The crossover behavior of the equation of motion (^) can be summarized in a F — L 
phase diagram which is shown in Fig. 9. On small scales L < the height fluctuations 
are bounded by the microscopic correlation length a±_ of the random forces. Here, the 
correlations (^ can be described by a smooth function ^{z). On larger length scales L^. < 
L < C, the system is in a critical state. The behavior of the interface is characterized by 
the critical exponents of the depinning transition. When crossing the Larkin length Lc, the 
second derivative A"(0) diverges under the RG flow and the fixed-point function A*{z) of 
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the correlator develops a cusp singularity at the origin. When going to length scales larger 



than ^ (for F > Fc), the behavior crosses over to the EW regime |jT9|. The full equation of 
motion (Q) has (for T = 0) an additional regime, namely the one of the KPZ equation . 
For D > 1 this regime is situated above the region of the EW equation. 

V. NUMERICAL RESULTS 

In this section we compare the analytical results for the critical exponents with extensive 
numerical simulations [0,|ll|]. Since our analytical values for the exponents (|61|) have been 



obtained by an expansion to first order in e, it is an open question whether or not the results 
(|61|) are applicable to e > 1. However, further work by Narayan and Fisher |^ suggests 
that the first-order results for the exponents ( and u are exact. We will see below that 
this is not consistent with our numerical simulations. The deviations between the analytical 
results and the simulations are much larger than the error bars on the numerical values but 
decrease with increasing D = 4 — e, thereby supporting our analytical calculation. 

A. Simulation of the continuum equation in D = 1 

We first simulate Eq. (0) in D = 1 with a discretization of the transverse coordinate only, 
X — > z, z{x, t) — > Zi{t). The random forces rii{zi) are chosen as follows: Each integer position 
Zi on a square lattice is assigned a random number t] between zero and one. For non-integer 
Zi the forces rii{zi) are obtained by linear interpolation [^. Finally, the ^-coordinates of 



rii{zi) in each column i are shifted by a random amount < Sj < 1, i.e., rii{zi) rji^Zi + Sj) 

At t = the interface is flat, Zi{t = {]) = Q. The interface configuration at t + At is 
calculated simultaneously for all i using the method of finite differences, 

Zi{t + At) = Zi{t) + At[F + Zi+i{t) + 2,_i(t) - 2zi{t) + mr]i{zi)\ . (62) 

Here m is a real parameter and periodic boundary conditions are used. 

We are interested in the critical exponents, so the threshold Fc has to be determined 
from a measurement of the velocity v{F) close to F^ where Eq. (|10D is valid. However, 
it turns out that the critical region is very narrow and that an accurate determination of 
v{F) inside that region is difficult, which seems to be the reason for contradictory results 



for 9 in the literature p8| , |30| , pT| . At the beginning of the interface motion the velocity de- 
creases and only saturates after a time when the roughness w ^ C,^ no longer increases. 
The time behaves as [F — Fc)~^'^ and therefore the simulations become increasingly 
time-consuming when approaching the threshold. In addition, the average velocity v is sys- 
tematically underestimated if the temporal fluctuations of the velocity become comparable 
with V itself. Thus, for F close to Fc very large systems have to be simulated. 

For m = 3 we estimated the threshold to be in the interval 1.41 < -F^ < 1.42 and for the 
velocity exponent we obtained an upper bound, 9{D = 1) < 0.45. 

To determine the roughening exponent (3 = (^/z, we measured the interface width 
w'^{L,t) = \[zi{t) — Zi{t)]'^) and the average position H(t) = {z{t)) in the transient regime. 
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t -C L^i for various forces in the interval 1.41 < F < 1.42. {L is the hnear system size, the 
overbar denotes the spatial average over i, and the angular brackets stands for the average 
over configurations of the disorder.) The best scaling of the width, w"^ ~ t^^, is achieved at 
F = Fc — 1.4135, with (3 = 0.88 ± 0.03, which is shown in Fig. 10. The average position 
H{t) increases with an effective exponent which is only slightly lower, (3h — 0.85. 

A possibility to determine both, jS and the roughness exponent ( is to consider the 
height-height correlation function C(r, t) = {[zi+rit) — Zi{t)Y). In the transient regime it 
scales as 

C(r,t) ~r2^^(rr^/"), (63) 

where \l/(|/) is a scaling function. Thus, the exponents ( and (3 can be obtained by a scaling 
plot of C(r, t) which is shown in Fig. 11. The best data collapse at F = Fc is achieved with 
({D = l) = zf3 = 1.25 lb 0.05 and = 1) = 0.88 ± 0.03. 



B. Simulation of the automaton model 

Since the simulations of the continuum equation (^) are computationally expensive, we 
also study a lattice model |jTO| of probabilistic cellular automata, which allows to determine 



effectively the critical exponents in higher dimensions as well. 

The automaton model is defined on a D + 1 dimensional hypercubic lattice where each 
cell [i, z] (with 1 < i < L^) is assigned a random force rji^z which takes the value 1 with 
probability p and rji^z = "~1 with probability 1—p. During the motion for a given time t the 
local force 

Mt)=Y.[^,{t)-Zi{t)]+mr^,,,^ (64) 

is determined for all i, where the sum is over nearest neighbors (j) only and m is an integer 
parameter. The interface configuration is then updated simultaneously for all i: 

Zi{t + 1) = Zi{t) + 1 if Vi>0 

Zi(t + 1) = Ziit) otherwise. (65) 

The parameter m measures the strength of the random force compared to the elastic force 
and the difference p — {l—p) = 2p — l determines the driving force. The growth rule specified 
by Eqs. (16^) and (|6^) can be derived from the continuum equation (^ by discretizing the 



transverse directions x, as well as time and the position of the interface z(x, t). In addition, 
a simple two-state random force rji^^ = ±1 is used which contributes to the computational 
advantage that one has to handle with integer variables only. 

The exponents ( and (3 are obtained in the same way as in Sec. V.A |1^. In addition. 



the static roughness exponent ( is determined by measuring the interface width of pinned 
interfaces for various system sizes L. The scaling w"^ ~ L^^ at threshold p = pc is shown in 
Fig. 12 {D =1,2 and 3). In D = 1 we obtain the same results for ( and (3 as for Eq. (|62D 
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but with smaller error bars. This supports the expectation that the continuum equation 
and the automaton model belong to the same universality class. 

The results for ( and f3 are shown in Table 1, together with the other critical exponents 
z = (/ (3, 6, and u. The correlation length exponent u is calculated from the scaling relation 
(|60|) which was shown to be exact (Details can be found in Refs. and fill)- The 
velocity exponent 6 is determined from the data of Fig. 13 (see Table 1). 



C. Self-organized automaton model (SOAM) 

The difficulty of the simulations of Eq. (|62|) and of the automaton model that the critical 
region is narrow can be circumvented by "self-organizing" the automaton model (|6^, (|65[). 
Then, the system is always at threshold without tuning a parameter This can be 



achieved by using an idea of Havlin et al. |^ and Sneppen p3|: Instead of updating the 



whole interface configuration simultaneously as in Eq. ( p5| ) one increments only that column 
j where the local force ( |6^ takes the maximum /max = niaxj[/j] among all fi |Tl| , 



Zi{t + 1) = Zi{t) + 1 if fi = /max 

Zi{t + 1) = Zi{t) otherwise. (66) 

To avoid ambiguities, the random forces 77,^2 are now chosen as real numbers in the interval 
< 1] < 1. A very similar model has been independently introduced by Roux and Hansen 
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We determine the roughness exponent ( of the SOAM by simulating interfaces according 
to the rules (0), (|66|) and measure the width during a sufficiently long time interval for 
different system sizes. The results ({D = 1) = 1.24 ±0.01, C(^ = 2) = 0.75 ±0.01, and 
C(-D = 3) = 0.34 ± 0.01 are consistent with the results of the automaton model. 



D. Comparison with analytical results 

We now compare the numerical results with the analytical exponents which are obtained 
by an extrapolation of the results (|6T|) of the e-expansion (see Table 1). It should be em- 
phasized that the equation of motion (|^) is not a useful description of an interface in a 
two-dimensional medium {D = 1) because higher gradient terms are relevant. This is in 
agreement with the unphysical result ({D = 1) ~ 1.25 > 1. Nonetheless, the considered 
interfaces are stable and self-affine because Eq. (^ and the automaton models (|64D , (|65|), 
and (|66D do not allow overhangs. We include the case D = 1 in Table 1 in order to compare 
the analytical and numerical results for Eq. (|^). 

The simulations support the analytical results: In D = 3 (e = 1) the analytical results 
deviate from the numerical values by only a few percent. In agreement with the expectation 
of an e-expansion in D = 4 — e the deviations become larger with decreasing D. The 
numerical results in Table 1 are consistent with the exact scaling relation (|TT]). However, 
the results for the roughness exponents are clearly inconsistent with the suggestion that 
( = [A — D)/3 is exact. Assuming that there are no higher-order corrections to C 0, the 
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numerical results indicate that there are non-perturbative corrections to (, a. possibility 
mentioned by Narayan and Fisher 0. 

The exactness of the result ( = {A — D)/3 would have meant that an interface at equi- 
librium [F = 0, T > 0) and a driven interface at threshold Fc have the same roughness 
exponent although the underlying physics is different. The simulations show that the driven 
interface is rougher than the interface in equilibrium. 

The Harris condition (|T3|) for a sharp depinning transition is fulfilled as an inequality by 
the numerical values. 



VI. SUMMARY 

In this paper we investigated the behavior of an elastic interface, which is driven by 
a homogeneous force F through a medium with short-range correlated random forces. It 
was made plausible that this problem can be described by the equation of motion (^) for a 
planar interface. It was found that = 4 is the upper critical dimension. For D > Dc and 
weak disorder, the interface is always smooth and there is no threshold for driven motion if 
the random force varies sufficiently slowly in space. However, we showed by considering a 
mean field theory that, a nonvanishing threshold Fc > can be produced by a random force 
with a cusp singularity (A'(O^) < 0) of its correlator A{z). Below Dc the interface becomes 
rough above the Larkin length L^. The correlation length ^ is a second important length 
scale, which diverges when approaching the critical point Fc. The results of the perturbation 
theory diverge with a power of ^. 

A simple three-parameter RG is not a way out of these divergences, because the second 
moment Q2 = A"(0) of the correlator and therefore the inverse mobility A diverge under 
renormalization aX L = Lc- Extending the space of renormalizable parameters to all mo- 
ments Q2n = A(^")(0) of the correlator, i.e. to the whole function ^{z), we found that the 
above divergences correspond to a cusp singularity of A{z) at the origin. With the interpre- 
tation Q2n = A*^^"'^(0"'") a consistent RG scheme could be constructed. The cusp singularity 
A'(O^) < produces a finite threshold Fc, similar to the mean field theory p3[. Above the 



threshold, the RG is stopped at the correlation length L = ^, above which one crosses over 
to a EW regime, where the random forces act independently on the moving interface. The 
behavior of the interface at different length scales L and driving forces F can be summarized 
in a phase diagram (see Fig. 9). 

The exponents in the critical regime Lc < L < ^ were calculated to first order in e = 4— D 
(see Eq. (pT|)). Although the result for the roughness exponent ^ ~ e/3 coincides with that of 
the equilibrium problem to first order in e, our numerical simulations show that the roughness 
at the threshold and zero temperature is larger than that in equilibrium (F = 0,T > 0). 
As expected, the velocity exponent 9 decreases when the dimension is decreased, i.e. for a 
rougher interface the onset of the motion at Fc is more abrupt. Our results for the correlation 
length exponent v fulfill a Harris-like criterion to observe a sharp transition as an equality 
for the first-order results and as an inequality for the numerical results. 
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APPENDIX A: VELOCITY 

In this Appendix the perturbative calculation of the velocity v of the interface is dis- 
cussed. We give two examples how typical integrals for the corrections to the velocity can 
be calculated. 



1. First-order correction 



The velocity v can be obtained by iterating Eq. (^21): The n — th order iteration for 
h(k, w) follows from (|22|) by using the {n — 1) — th order iteration on its rhs. (The zeroth 
order iteration of h vanishes.) After truncating this procedure at the order which is needed, 
the resulting equation is averaged over the disorder, which gives on the Ihs {hCk,^)) = 0, 
while the rhs is a function of v. For example, the lowest-order correction (Fig. 4) to 
Vq = F/X is obtained from the second-order iteration 

F-A„., Iff A, A„(P) _ 



\v \v Jq Jp — iXvq 

where we have included the Fourier transform A||(p) of the correlator along the interface 
(see Eq. (g|)). 

First, we make Ji real, perform the angular integration and take the correlator A||(p) 
into account by truncating the p- integration at l/ay: 

where K]^^ = 2^^^tx^/'^V{D /2). For D > A this integral converges for all v as long as ay is 
finite. For D < 4, however, the limit f — > causes a divergence at small p. If we set a\\ = 0, 
the integral (A2) diverges at large p when 4 — Z) = e ^ 0. 



After introducing the new variable y = p/ JXvq/'y, we get 



i,^p,„r"i/'''^,l'^""^yf^,. (A3, 

The divergence with f — is now outside the integral and we can set the upper limit to 
infinity and thereby getting the desired leading term proportional to 1/e: 



h ^ A"(0) / dy y^-\ (A4) 

7^ Ji 
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where Eq. ( p6D was used. Hence one gets the result 

^ cQ.hl = 5A(iVA, (A5) 

which is the same as Eq. (|37|). 

However, for very large v , a\\, i.e. yXvqJ^ can no longer be regarded as a lower 

cutoff for the p-integration. The latter has to be truncated at an upper cutoff l/ay. In this 
case one can introduce the new variable q' = Xvq/'jp'^. Then the integrals in Eq. (|A2|) can 
be carried out: 

- " ^N'T^-A(0)T^4r- (A6) 



As expected, the correction to F — Aw vanishes for large f , where Eq. is of EW-type. 

2. Second-order correction 

Next we extend the calculation of the velocity to second order. To this end Eq. ( |22D 
is iterated as described above. The second order diagrams for the velocity are shown in 
Fig. 14. Since we want to use the perturbation theory to construct a renormalization group 
expansion in D = 4 — e we only calculate the leading order terms in e. The expression for 
diagram 14a) is given by 

I I I i A, A,, A||(p)A||(p') t'q'q' 



q Jq' Jp Jp' 



X 1 1 1 (AT) 

7P^ + iXvq 7p'^ + iXvq' 7(p + p')^ + iXv{q + q') 

The integrals will be evaluated below. The result is 



The calculation of diagram 14b) is very similar and gives the same contribution as ( |A8D . 
Diagrams 14c) and 14d) correspond to a renormalization of the response function and also 
contribute the term ([A8|). The sum of the diagrams 14e)-g) is equal to zero, diagram 14h) 
cancels with 14j) and diagram 14i) cancels with 14k). If we collect all diagrams up to second 
order we finally obtain 

2 



Ae// - A 



, 5A« ,AA«\ 



(A9) 



As an example for the evaluation of second-order diagrams we calculate the contribution 
of diagram 14a). The integrals over q and q' do not vanish only if the powers of q and 
q' are even. Thus I2 is 



/ / A, A,,A,|(p)A||(p')g^g'^t; 
jp jp' 



q Jq' Jp Jp' 
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X 



pV' + p^(p + p')' 



[72p4 + {Xvqy] [72p'4 + (Xvq>)2] [^2(p + p/)4 + ^Xv{q + q')^ 



(AlO) 



A, A,, q' q" V h. (All) 



q Jq' 



We consider here only the first term with p^p'^ in Eq. ( |A1CI| ) because the second one cancels 
with a term of the diagram in Fig. 14b). We have seen that the integration over p gives 
the important leading divergences, while the g-integration does not yield divergent factors. 



Thus, for simplicity we set in Eq. ( |A11| ) q = q', and introduce new variables y = p/ JXvq/'y, 



and y' = p' / y )^vq/'~f and treat the limits as in the calculation of the first order correction. 
Then we obtain 



(27r)^A (27r)^ y^ y'^ (y + y') 



(A12) 



where the range of integration refers to the absolute values of y and y'. Again, the important 
divergence is at large y and y' , so that we can introduce a smooth cutoff for small y and y' 
and first perform the ^/'-integration and set the upper limit to infinity: 



y (y= + 1) 



k 



with 



K'{y" + i)[{y + y'f + ir ^^^^^ 



This D-dimensional integral can be evaluated by introducing a Feynman-parameter x ||36 



4 = j^^ r(3) dx X {(1 - x){y'^ + 1) + a;[(y + y')' + 1]}"' • (AM) 

Using the general formula 

f\ 2^9 2i-a Kp r (f ) F (g - f ) ,21(^/2)- 

/ [y + 2yy + m ] = — ' .\ m -y , (A15) 

one obtains 

/' = :^r(D/2)F(3 - D/2) [' dx X \y\x - x^) + 1]^''^^^''' . (A16) 
2 Jo ^ J 

The y-integration is also performed by introducing a Feynman-parameter x'. 

r\D/2)r{e) 



~ V 2~ 

X dx [' dx' X f^n^^i _ _ 1] + i|(^/2)-2 _ ^^^^N 
Jo Jo 
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With e/2 = e' — ^ this integral is logarithmically divergent. The pref actor of this logarithmic 
divergence is the same as the prefactor of the 1/e term for small e. Thus, we consider only 
the integral Eq. ( [A17| ) for e' = e/2 = and calculate the desired prefactor. The integrals are 



then elementary and the prefactor is calculated to one. Then we obtain with r(2e') ~ l/(2e') 

(A18) 



With Eq. PTT[ ) the result (^) follows. 



APPENDIX B: VERTEX CORRECTION 

In this Appendix we outline the calculation which leads to the renormalization of mo- 
ments resulting in Eq. (|6). The diagrammatic representation of the bare moments Q2n 



is shown in Fig. 7a) , where all incoming and outgoing momenta conjugated to x are set 
to zero. According to Eq. (^) the analytical expression associated with the diagram a) 
is exactly Q2n- The diagrams contributing to the renormalization of the moments up to 
the one-loop order are shown in Fig. 7b)-d). In each of these diagrams, one has only one 
integration over the internal momentum conjugated to the space argument of the interface 
height /i(x). This integration has the form /pG'q(p) and yields for all diagrams under con- 
sideration the same factor c^^e. Additionally to this integration there is also an integration 
over the momenta qi and q2 conjugated to the argument of the disorder correlator A(z). 
The momenta qi and q2 are associated with the dotted lines in one-loop diagrams in Fig. 7. 

Let us consider the diagram b). The number of arrows starting at crossing points of 
dotted lines with the continuous line is n for both upper and the down continuous lines. The 
distribution of arrows between the both crossing point along a continuous line is arbitrary. 
The sum of the arrows should obey the condition rii + n2 = n. According to Eq. (|2^) the 
factor (±ig)''" is associated with each crossing point with m being the number of lines starting 
at the point under consideration. The sign -|- has to be taken if the dotted line belonging to 
the crossing point goes out from the point, — if the dotted line goes into the crossing point. 
In the following we designate —q by q. Summing over all possible distributions of arrows 
between the two crossing points of the top line in the diagram b) we get the factor 

i''^\i{qi + q2r. (Bl) 
The factor associated with the distribution of arrows on the down line in b) is 

^"+'gl(gl + g2r, (B2) 

so that the factor associated with the whole diagram b) is 

(-ir+Vgi(gi+g2r(gi+g2r. (bs) 

Analogously the factor associated with the diagram c) is given by 

{-lT^'qMqi + q2T{qi + q2T- (B4) 
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The factor associated with b) and c) is the sum of (p3|) and ( p^) 

2n+l 
k=l 

where Cj = (l) are the binomial coefficients. The diagram d) in Fig. 7 gives the factor 

- qr'-'ql (B6) 

which compensates the last term in Eq. (|B5|). Thus, the part of the diagrams b) - d) 
depending on qi and q2 is given by 

2n+l 

E +igf (B7) 

k=l 

To obtain the analytical expression associated with b)-d) we have to multiply ( p7|) with 
A(gi)A(g2) and integrate over qi and g2- Taking into account that only the even moments 
of the disorder correlator A{z) are non zero we get from (E?) 

n 

E ^2fc^ Q2kQ2n-2k+2- (B8) 
k=l 

The whole correction of b)-d) is obtained from ( P8| ) by multiplying the latter with the factor 
c^o/e. Finally, the contribution from the diagrams a)-d) in Fig. 7 is 

n 

Q2n + E ^2k^^Q2kQ2n-2k+2C^Q/ ^- (B9) 
k=l 

The latter gives the renormalization of the moments of the disorder correlator up to the 
one-loop order. 



APPENDIX C: FUNCTIONAL FLOW EQUATION 

In this Appendix it is shown how the vertex corrections of the moments Q2n in Eq- (^61) 
can be put together in a functional flow equation for A{z). 

The flow equation can be formally obtained by differentiating Eq. (|26| ) 

d A{z) _ ^ jtz)^- dQ^n 

d InL (2n)! dL ' ^ ^ 

where Q2n is identified with the effective moments Q2n,eff which depend on the upper cutoff 
L = ^0- Differentiating Eq. (^61) and writing Q2n,eff = Q2n one gets 

L = C-^' E C!2j~^^Q2jQ2{n-j+l)- (C2) 

To construct a functional fiow equation the rhs of Eq. ( |CI| ) has to be expressed as a 
function of A(z). To this end one first extends the summation over j to < j < n and 
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subtracts the term with j = separately. Introducing the new index k = n — j and using 



/2j — C2(n-j)+i' ^he rhs of Eq. (|C1| ) can be written as 



n=0 ^0 



H C'2j+1 Q2i+2Q2n-2j " <5o /r,^N, Q2n+2, 

n=0 l^'^J- 



(C3) 



where j is written for A; and the factor cL^ is skipped. With Eq. ( p6D the second term can 
be expressed as 



92 



d{iz) 



rA(^)A(0). 



(C4) 



The first term in (|C3|) can be written as 
d 



(2n-2j)!^"' 

z 

Writing down the sums over n and j term by term, one easily sees that 

d 



ditz) 1 ^i^z') fr{ (2j)! 



00 t-.2n 



,_,n?o (2-)! 



Q 



2n 



(C5) 
is equal to 

(C6) 



The range of the summation over j can be extended to the vanishing term j = 0. Using Eq. 
J6D can be identified with 



d 



d 



A(z) 



1^ 
2d^ 



A\z). 



(C7) 



d{iz) [ d{iz') 

Now we can put the first and second term of Eq. ([C^ ) together and insert it in Eq. ( |Cl| ) : 



d A{z) _ 

d InL dz"^ 



-A\z) - A{z)A{0) 



(C8) 



which is the functional flow equation 
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FIGURES 



FIG. 1. Interface veloeity f as a function of the driving force F. 

FIG. 2. The random force correlator A{z) for (a) random-fields, (b) random bonds, and (c) 
charge-density waves. 

FIG. 3. The diagrammatic representation of Eq. (22) (see text). 

FIG. 4. First-order correction to the velocity (see Eq. (29)). The average over the quenched 
randomness is represented by the connection of two outgoing dotted lines, which build up a "loop" . 

FIG. 5. The response function G(k, w), represented as a double line with an arrow, and its 
one-loop expansion (see Eq. (33)). 

FIG. 6. Lowest-order diagram to the interface width (Eq. (39)). 

FIG. 7. Diagrams illustrating the renormalization of the moments of A(2:). The diagrams 
shown contribute to the renormalization of Q2n- In diagram b) and c) the number of outgoing 
lines have to fulfill the condition m + n2 = + = n. 

FIG. 8. When starting with a smeared delta function (full line in figure on the Ihs) the second 
moment Q2 = A"(0) is negative. Under the renormalization on scales L < Lc the curvature at 
the origin becomes greater and A{z) becomes more peaked (dotted line). At L ~ Lc the curvature 
diverges and the fixed-point solution has a cusp singularity (figure on the rhs). 

FIG. 9. F — L phasedigram (see text). 

FIG. 10. Scaling of the average position H{t) (triangles) and interface width 'w{t) (circles) at 
threshold Fc from a simulation of Eq. (62) in D = 1. The system size is L = 16384, At = 0.125, 
and the data are averaged over 30 realizations of the disorder. 

FIG. 11. Scaling plot of the height-height correlation function C{r,t) from the same runs as 
in Fig. 10. The same plotting symbol is used for data at a given time t. 

FIG. 12. The width of pinned interfaces as a function of the system size L from 
simulations of the automaton model for D = 1, 5 = 1, p = — 0.8004 (open circles), 
D = 2, g = 6, p = pc — 0.74446 (filled circles), and D = 3, g = 8, p = pc^ 0.63165 (squares). 
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FIG. 13. Interface velocity f as a function of the driving force p — Pc for D = 1, g = I, 
L < 262144 (open circles), D = 2, g = 6, L"^ = 1024^ (filled circles), and D = 3, g = 8, L < 110^ 
(squares) . 

FIG. 14. Second-order corrections to the velocity. 
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TABLES 



iLixponent 


U 


1 J. " 1 

analytical 


numerical 




-1 

1 


1/3 


0.25 zb 0.03 


velocity 


2 


2/3 


0.64 ± 0.02 




i 


U.ODO.. 


0.84 ± 0.02 




1 


1 


1.25 lb 0.01 


roughness C, 


2 


2/3 


0.75 ± 0.02 




3 


1/3 


0.35 zb 0.01 




1 


3/4 


0.88 zb 0.02 


roughening /? 


2 


0.4286.. 


0.48 zb 0.01 




3 


0.1875 


0.20 zb 0.01 




1 


4/3 


1.42 zb 0.04 


dynamical z = C,/ (5 


2 


1.555.. 


1.56 zb 0.06 




3 


1.777.. 


1.75 ± 0.15 




1 


1 


1.33 ± 0.02 


correlation length u 


2 


3/4 


0.80 zb 0.01 


u = l/{2-0 


3 


3/5 


0.606 zb 0.004 



TABLE I. Analytical and numerical results for the critical exponents for interface dimensions 
D = 1,2 and 3. The analytical values are obtained by an extrapolation of the results (61) of the 
e-expansion. 
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